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We analyze the effects of the electronic coupling to bosonic modes in a d-wave superconductor. 
The role of the scattering due to boson on the momentum transfer between electronic states in the 
Brilloine zone is addressed. We consider specific examples of Big phonon, breathing mode phonon 
and spin resonance at {•K,Tr). The Fourier spectrum of the energy derivative local density of states 
(LDOS) is calculated. To properly calibrate the effects of different modes we fix the quasipartilce 
renormalization at specific momentum points. It is found that the Big mode with highly anisotropic 
momentum-dependent coupling matrix element gives rise to well definded features in the Fourier 
spectrum, at the energy of mode plus gap, with a momentum transfer along the Cu-0 bond direction 
of cuprates. This result is in a striking contrast to the cases of the coupling to other modes and also 
to the case of no mode coupling. The origin of this difference is explored in detail. A comparison 
with the recent STM experiments is briefly discussed. 

PACS numbers: 74.25.Jb, 74.50.+r, 74.20.-z, 73.20.Hb 



I. INTRODUCTION 

Determining the nature of single particle excitations 
is of fundmental importance in our understanding of 
the superconductivity in high-Tc cuprates. To address 
this issue, a number of spectroscopies have been exten- 
sively used, including the angle resolved photoemission 
spectroscopy (ARPES) and tunneling. The salient fea- 
tures observed in ARPES include: (a) Near the (tt, 0) 
(M) point in the Brillouin zone, the spectral func- 
tion in the superconducting state shows an anomalous 
line shape, the so-called "peak-dip-hump" structure>i*2i^ 
(b) Near the d-wave node of the superconducting gap, 
the dispersion shows a characteristic "kink" near 50-70 
j^gYj4j5j6j7j8 i^gcent ARPES experiments with improved 
resolutior&iSiiiii^ have revealed another "kink" in dis- 
persion of the antinodal electronic states, near the M 
point. An unusual spectral dip- hump features similar 
to the ARPES spectrum have also been observed in the 
tunnehng datai^^i^'^i^^i^''i^''i^^i^^i^° All these features were 
suggested to indicate that the electron self-energy renor- 
malization could be due to the electronic coupling to 
a bosonic mode. Two main scenarios have been pre- 
sented to explain the experimental data. On one hand, 
the antinodal renormalization is found to be strongly en- 
hanced below Tc AiSiii Such a strong temperature depen- 
dence and the dominance of the coupling strength near 
the M point can be thought of as evidence for the cou- 
pling of electrons to the 41 meV spin resonance mode 



(of electronic origin) 
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As seen by in- 



elastic neutron scattering experiments in most of the 

CUprateS,3"i31i32,33,34,35,36,37,38,39,40,41,42^43,44,45 ^^^^ 

mode intensity substantially turns on below Tc (even 
though some intensity might be present in a normal state) 
and has a well-defined momentum of (tTjTt). This sce- 
nario of the electronic coupling to spin resonance mode 
has also been used to explain the tunneling spectra in 
planar tunnel junctions. On the other hand, it has been 
suggested that a significant electronic coupling to the 
half-breathing in-plane Cu-0 bond stretching phonon or 
to the out-of-plane out-of-phase O buckling Big phonon, 
with an energy of approximately 70 and 35 meV, respec- 
tively, might be responsible for the dispersion anomalies 
at the nodal^ and antinodal directions respectively. 
These two phonon modes have shown strong lineshape 
renormalizations with doping and temperature in Ra- 
man and neutron measurements 4Si^£i^242iS24^iiS4S The 
advantage of this scenario is that it could naturally ex- 
plain the band renormalization effect in materials where 
no spin resonance mode has been detected, in the normal 
state, and in the deeply overdoped region where the spin 
mode is neither expected nor observed. To be consistent 
with ARPES data, this scenario requires the electron- 
phonon interaction to be highly anisotropio^^i^ and its 
impact on the electrons to be strongly enhanced in the 
superconducting state. 

The nature of the involved bosonic modes, being 
phononic or electronic and their role in the mechanism 
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of superconductivity remains controversial. Detailed en- 
ergy and momentum spectoscopy of the relevant bosonic 
modes might be very helpful in understanding the mech- 
anism of superconductivity in high- Tc cuprates. One 
of the direct spectroscopies that allows energy deter- 
mination is Inelastic Electron Tunneling Spectroscopy 
(lETS). It is a well-established and powerful tool that 
allows the measurement of the characteristic energies of 
extended modes. Examples of applications of this tech- 
nique, among many, include measurements of molecu- 
lar stretching and vibrational modes in metal-insulator- 
metal tunnel junctionsj^SiS observation of the collective 
magnetic resonance in tunneling in the superconducting 
state of high- Tc materials^fliSi and observation of the tun- 
neling features at energies that correspond to the phonon 
peaks, as seen in planar tunneling into superconduc- 
torsiS 

The lETS directly measures excitation energies. When 
electrons scatter off a collective mode, a contribution 
to the electron self energy occurs above a correspond- 
ing threshold value of the frequency determined by the 
mode frequency. Thus, in a tunneling experiment, for 
bias voltages exceeding the threshold, electrons can ex- 
cite the mode. This additional scattering channel leads 
to a step in the density of states (DOS) and in the tun- 
neling conductance. Low temperatures are required to 
avoid thermal smearing of the step in the conductance. 
The crucial quantity that reveals inelastic peaks is a sec- 
ond derivative of tunneling conductance with respect to 

,2 T 

bias voltage: -^^^ieV). Peaks in this quantity are shown 
to be connected to the energies of the modes, e.g peaks 
in phonon Density of States 

In high- Tc materials the energies of a number of col- 
lective modes, like phonon and spin modes are close. For 
example the Big phonon mode has typical energy of 36- 
40 meV and spin resonance mode has an energy that 
ranges between 35-40 meV depending on doping. Hence 
discrimination between different modes based on only the 
energy of the observed mode is a challenge. Aside from 
energy resolved features, it would be useful to come up 
with the measurement that would allow one to measure 
typical momenta involved in the electron scattering in 
cuprates. The ARPES is one such spectroscopy. 

Scanning Tunneling Spectroscopy (STM) is another 
technique that would allow one to resolve the momentum 
transfer between different electronic states. Recently, the 
Fourier-transformed scanning tunneling microscopy (FT- 
STM)SiS£ was introduced to map out the Fermi surface 
and the momentum dependence of the d-wave energy gap. 
The results are in good agreement with the ARPES data. 
While the ARPES provides the information about the 
energy dispersion of single-particle electronic states, the 
FT-STM probes the scattering processes between states 
with different in plane momenta. 

In this paper, we show how another technique, the 
Fourier Transform Inelastic Electron Tunneling Spec- 
troscopy (FT-IETS), would allow simultaneous momen- 
tum and energy resolutions of the tunneling electrons. 



As such, the technique might be useful to address the 
role of different modes in cuprates. Motivated by the 
progress of elastic FT-STM technique, earlier few of us 
(JXZjQS, AVE) have suggested an FT-IETS technique, 
Ref. |63, in a specific model that describes the Fourier 
transformed local density of states (LDOS) in a d-wave 
superconductor with the electronic coupling to the (tt, tt) 
spin resonance mode. Here we go beyond the previous 
analysis and address the question of what consequences 
the electron self-energy renormalization will have on the 
tunneling characteristics for a number of collective modes 
broadly considered in the literature: the Big and breath- 
ing phonon modes, as well as the (tt, tt) spin resonance 
mode. The central quantity we will focus on in our anal- 
ysis is energy derivative of the FT LDOS, p'{q,E). This 
quantity corresponds to (q, e V) measured by the FT- 
STM, where / is the local tunneling current and V the 
voltage bias. Here we are not concerned with the mecha- 
nism of the superconductivity in the cuprates. Instead we 
assume from the outset a d-wave channel effective pair- 
ing interaction and study the additional effects due to the 
electronic coupling to various bosonic modes, including 
the Big and breathing phonon modes, and the (tt, tt) spin 
resonance mode. A comparison of the calculated momen- 
tum transfer structure with the FT-STM measurement 
may shed new light on which type of bosonic mode the 
electron excitations are coupled to strongly. 

The rest of the paper is outlined as follows: In Sec. ITU 
we develop a theoretical model in which the electrons 
are coupled to bosonic modes. We consider three types 
of modes: the Big and breathing phonon modes with 
possible coupling matrix elements, and the (tt, tt) spin 
resonance mode. Additional weak disorder is used as 
a marker so that the momentum transfer can be inves- 
tigated. In Sec. mil the numerical results for the FT 
spectrum of the energy derivative LDOS arc presented. 
Sec. IIVI contains some concluding remarks. 

II. THEORETICAL MODEL 

We start with a BCS-type model to describe two- 
dimensional electrons with a d-wavc pairing symmetry, 
which is relevant to high-temperature cuprates: 

T-Cbcs = Xl^k4^ck^ + ^(Akcj^^cL^^ 4- A;;c_kiCkT) , 

k,(T k 

(1) 

where cj^^ (cko-) creates (annihilates) a conduction elec- 
tron of spin (T and wavevector k. The quantity ^k is 
the normal state energy dispersion. We adopt a six- 
parameter fit to the band structure used previously for 
optimally doped Bi-2212 systemsjS having the form 

^k = — 2ti{cos kx + cos ky) ~ 4t2 cos kx cos ky 
— 2i3(cos 2fc^ + cos 2ky) 
—4^4 (cos 2kx cos ky + cos k^ cos 2ky) 
—4:t5 cos 2kx cos 2ky — /i , (2) 
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where ti = 1, = -0.2749, h = 0.0872, = 0.0938, 
is = —0.0857, and ^ = —0.8772. Unless specified explic- 
itly, the energy is measured in units of ti hereafter. 

As in previous works ^SSi^iiSiSi we assume that an ef- 
fective d-wave pairing interaction has pre-existed from 
certain many-body effects. Therefore, even before the 
electronic coupling to the bosonic excitations, the d-wave 
superconducting order has already been established and 
the corresponding order parameter is given by: 

Ak = ^(cosfcj, - cosfcy) . (3) 

The additional renormalization comes from the addi- 
tional electron-bosonic mode interaction. 

By introducing a two-component spinor operator, one 
defines a matrix Green's function in the Nambu space. 
The bare Green's function in the assumption of a real 
d-wave pair potential reads 

Here w„ = (2n + 1)tiT is the fcrmionic Matsubra fre- 



quency. 



A. Coupling to the collective modes 

As mentioned before, there are many bosonic collective 
modes existing in cuprates. Here we focus on the out-of- 
plane out-of- phase buckling {Big) and the in-plane half- 
breathing phonon modes associated with the motion of 
the O ions, and the (tt, tt) spin resonance mode. 

We model the electronic coupling to the phonon modes 
by 

Hei-ph = — ^ 5Z q)Ck-Hq,<TCk^^>^,q ' (5) 

* ^ k,q 

where Nl is the number of lattice sites, Al,_c^ = bu<^+b\, _q 
with fctq (fei^q) creating (annihilating) one phonon of type 
V (representing Big or breathing mode) and wavcvector 
q. We consider two types of the coupling matrix ele- 
ments. The first type are both q and k dependent: 



5i3ig(k,q) 



.90 



{0a;(k)(/)2;(k -I- q) cos{qy/2) - (f)y(k)(f)y(k + q) cos{qx/2)} 



.(k,q) = .go ^ {(f)b{k + q)(j)a{k)cos[{ka+qa)/2]~ (pbik)(l>a{^ + ci)cos{ka/2)} 

a—x^y 

I 



(6) 
(7) 



where M(q) = [cos'^{qx/2) + cos'^ {qy/ 2)] /2, and 



VVk 




(8) 
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Ck^y,k ixy,\ztxM 


(9) 


1 
A^ 


f2 _ 4-2 

Sk ''xy.k. 1 


(10) 



with AAk = {[^l - tlyj"^ + [Ckt^.k - txy,kty,k? + 

[£.kty,w - txyMtx,k]'^Y^^, ta,w = -2ii sin(fc„/2) and 
txy.k ~ — 4^2 sin(fc2:/2) sin(fcy/2). The k dependence 
has been argued to be crucial in the interpretation of 
ARPES dataii^iS The second type has only q depcn- 

,21.54.63.64.65.66 



|gB,Jk,q)p = |5oP 
|gb.(k,q)p = |go|2 



2 1x 



2 My 



■ 2 I'-i 

sm ' 



cos 



■ 2 / Qv 
sm ' 



(11) 
(12) 



This type of q has been the focus before in the context 
of a dx2-y2-weLve pairing mechanism. For convenience of 



notation, we refer to the phonon modes with the first 
type of coupling as Big-l and br-1 ones, while those with 
the second type of coupling as Big-ll and 6r-H ones. 

The electronic coupling to the (tt, tt) spin resonance 
mode is modeled as: 



90 
2Nt 



(13) 



where S is the spin operator for the (tt, tt) mode, cr is the 
Pauli matrix vector. 

We calculate the electronic self-energy due to the 
electron-bosonic excitation coupling up to the second or- 
der in the coupling matrix elements. For the electron- 
phonon coupling, the self-energy is given as: 



E(k;iw„) 



T 



^^9^0^- (I, q)5i^ (k, -q) 



xX'^(q;if7,„)f3^o(k - q; - i^m)T3 , 

(14) 

where flm = 2rmTT is the bosonic Matsubra frequency, 
is the third component of the Pauli matrix in the 
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Nambu space, the quantity X',y(q;in,„) is the Fourier 
transform of the phonon Green's function I?j,(q;T) = 
— {Tr[A^,q{T)A^,^q{0)]) and is taken as 



1 



(15) 



with Qi, being the resonance frequency of the phonon 
modes. 

For the electron- (tt, 7r)-resonance spin-fluctuation cou- 
pling, the self-energy is given asiSi 



I](k; iujn 



Ml 

4Nl 



q 0„ 



(16) 

where x(q; ^^m) is the Fourier transform of the spin- 
spin correlation function, x(q; r) = — (7V[<5'q(T)S'iq(0)]) 
dynamical spin susceptibility. We treat the susceptibility 
in a phenomenological formi^ 



X(q; iilrn) = - 



/(q) 



1 



1 



(17) 



Here the spin resonance mode energy is also denoted by 
Hq. The quantity /(q) describes the momentum depen- 
dence of the mode and is assumed to be enhanced at 
the Q — (tt, tt) point. Using the correlation length ^sf 
(chosen to be 2 here), it can be written as 



/(q) = 



1 



(18) 



The form of the dynamic susceptibility as given by 
Eq. (|17|l is especially suitable for the optimally doped 
YBa2Cu306+i/ (YBCO) compounds in the superconduct- 
ing phase, where the observed neutron resonance peak is 
almost resolution-limited in energy and fairly sharp in 
wavevector. The resonance peak in BSCCO is broad- 
ened in both energy and wavevector. In addition, given 
that the peak in BSCCO is still quite sharp in energy, 
we expect that the energy broadening of the resonance 
mode is not important for the present study. We have 
also neglected the incommensurate peaks seen in the 
inelastic neutron scattering experiments in YBCO (the 
part that disperses "downward" away from the resonance 
peak) (MiSLSSiSSiZS since their spectral weight is signifi- 
cantly smaller than that of the resonance mode. 

The dressed electron Green's function due to the 
renormalization effect of bosonic excitations, is given by: 



Q ^ (k; iuJn) = Ga ^ (k; iuJn) - S(k; iu;„) 



(19) 



To study the momentum transfer between the bosonic 
excitation renormalized electronic states, additional im- 
purities or defects are required to scatter the electrons. 
The scattering from impurities is described by, 



(20) 



where Ui is the strength of the zero-ranged impurity scat- 
tering potential at the l-th site. For simplicity, we con- 
sider only the case of nonmagnetic scattering and assume 
the scattering potential from all these impurities are iden- 
tical, i.e., Ui = Uq. The full Green's function satisfies the 
following equation of motion: 

G{i,j;iujn) = g{i, j;iuj„)+^Uig{i,l;iLUn)f3G{l, j;iujn) ■ 

(21) 

Since we are most interested in effects of the electron col- 
lective mode coupling, it is desirable to keep as small as 
possible the disturbance (i.e., quantum interference, for- 
mation of virtual resonance etc.) arising from the quasi- 
particle scattering off the impurities themselves. This 
suggests to consider a dilute concentration of weak impu- 
rities. In this limit, the Born approximation is applicable, 
and one arrives at: 

G{i,j; iLOn) = G{i,j; iujn)+Uo ^^G{i, I; iujn)T3G{l,j; iwn) ■ 

Due to the impurity scattering, the correction to the 
local density of states (LDOS) at the i-th site, summed 
over two spin components, is 



2Uo 



(23) 



B. Fourier Transform 

The local density of states is proportional to the local 
differential tunneling conductance (i.e., dl/dV). To look 
into the renormalization effect of collective bosonic exci- 
tations in the STM, the energy derivative of the LDOS, 
corresponding to the derivative of the local differential 
tunneling conductance (i.e., cfi I /dV'^), is more favorable 
to enhance the signal. For a fixed value of energy, one 
first gets a set of Sp'{i, E) (the prime / means the energy 
derivative) in real space, and then performs the Fourier 
transform: 

5p'{ci,E)^Y.^p'{v,,E)e-'^-^^ , (24) 

i 

to obtain a map of the Fourier spectrum in q space, 



(25) 



P{ci,E) = \5p'{ci,E)\ 

One can also prove that the relation between (5p(q, E) 
and that due to a single impurity 5psingie{(\^ E): 

5p'{ci,E)=F{ci)5p',„,giM,E) , (26) 

where i^(q) is the form factor for the spatial distribution 
of weak impurities [F{c\) = 1 for the case of a single 
impurity] . 
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Big-1 mode br-1 mode -B19-II mode br-11 mode (tt, tt) mode 
1.05 1.48 0.75 0.75 2.30 



TABLE I: The calibrated value ol the coupling strength go 
for different type of collective modes. 

III. RESULTS AND DISCUSSIONS 

For the numerical calculation, we take the supercon- 
ducting energy gap Aq = 0.1, the frequency of all col- 
lective modes fio = 0.15. A quasiparticle lifetime broad- 
ening of 7 = 0.005 is used. A weak impurity scatter- 
ing strength Uq = 0.1. We take a large system size 
Nl = 1024 X 1024 to achieve the high-momentum and 
energy resolution. The Fourier spectrum P(q, E) is then 
constructed from 5p'(ri; E) within a given window of size 
61 X 61 centering the single impurity. We choose the cou- 
pling strength for all three types of collective modes in 
such a way that at the Fermi energy £' = 0, the frequency 
renormalization factor Z appearing in the self energy 

S(k;icj„) = iw„[l-Z(k;iw„)]fo-fx(k; jcj„)f3-H$(k;icj„)fi , 

(27) 

has the same real-part value for the Sig-I, 6?'- 
II phonon modes, and (tt, 7r)-rcsonance spin fluctuation 
modes at the M point, while for the hr-1 phonon mode 
at the wave- vector (7r/2,7r/2). The calibrated value of 
the coupling strength go for all these collective modes is 
summarized in Table 




FIG. 1: Density of states (left column) and its energy deriva- 
tive (right column) as a function of energy for a clean d-wave 
superconductor with the electronic coupling to Big-I, br-I, 
and (tt, tt) spin resonance modes. The case of no mode cou- 
pling is also shown for comparison. 

We present in Fig. ^ the results of the DOS and its 
energy derivative as a function of energy for a clean (i.e., 
Uq = 0) d-wave superconductor with the electronic cou- 
pling to the Big-1, br-I, (tt, tt) spin resonance modes. For 
comparison, the DOS for the case of no mode coupling 



is also shown. When there is no electron-mode coupling, 
there is a van Hove singularity peak appearing outside 
the superconducting gap edge. When the electrons are 
coupled to the Big and (tt, 7r)-spin resonance modes, the 
van Hove singularity peak is strongly suppressed. In- 
stead, one sees a dip structure following the coherent 
peak at the gap edge. The distance between this dip 
and the coherent peak defines the resonance energy f2o- 
However, there is very little suppression when the elec- 
trons are coupled to the br-1 mode. The planar tunnel- 
ing experiments indeed observed the peak-dip structure 
rather than the peak-peak (van Hove singularity) struc- 
ture. For the band-structure parametrization we have 
adopted, this implies that the electronic coupling to the 
collective modes must exist, and the Big and (7r,7r) spin 
resonance mode are the most promising candidates. Un- 
fortunately, as shown in Fig. ^ the dip structure due to 
the coupling to the Big and (tt, tt) spin resonance mode 
is almost identical. It would be very challenging to dis- 
tinguish between these two modes in the planar tunnel- 
ing experiments, which is measuring the momentum av- 
eraged DOS. Therefore, we propose to look at the mo- 
mentum transfer structure between the band renormal- 
ized states, which can be measured by the FT-STM. To 
achieve this goal, we need to have a signal strong enough 
to be detectable in STM experiments. The derivative of 
the DOS, p'{E), would serve the purpose. As shown in 
the right column of Fig. ^ when the electrons are coupled 
to the Big and spin resonance modes, there is a strong 
peak structure a,t E — — (Ag+ilo) exhibiting in the p'{E) 
spectrum, which has a one-to-one correspondence to the 
dip structure in p{E) itself. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



Br-M (it.it) No mode 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



FIG. 2: The Fourier spectral weight of the energy derivative 
of the LDOS at E — — (Ao + Slo) for a d-wave superconductor 
with the electronic coupling to the Big, br-I, _Big-II, br-II, 
spin resonace modes. For comparison, the quantity is also 
shown for the case of no mode coupling. 

In Fig. 13 we present the results of the Fourier spec- 
trum of the derivative of the LDOS, p'{q,E), at the en- 
ergy E = — (Aq -I- ^q) for a d-wave superconductor with 
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the electronic coupling to the collective modes. For com- 
parison, the same quantity is also shown (last panel) for 
the case of no mode coupling. Note that the case without 
the mode coupling, the energy Qq has no special meaning 
in the context of the electronic properties, and the energy 
E = — (Aq + fio) is chosen merely for comparison to the 
case of mode coupling. The main results are: For all 
cases, there are strong intensity at the large momentum 
transfer near (tt, tt). For the cases of the br-l coupling 
and no mode coupling^ similar feature at very small 
momentum transfer is obtained, relating to the fact that 
the DOS spectrum in two cases (see Fig. ^ is similar to 
each other. For the cases of the coupling to the -Big-II 
and br-11 modes, there are intensity peaks at a finite mo- 
mentum transfer along the diagonals. For the case of the 
coupling to the spin resonance mode, no strong feature is 
obtained at the intermediate value of momentum trans- 
fer, which is consistent with our previous calculation for 
this specific casei^ For the case of the coupling to the 
Big-1 phonon mode, there exists not only the intensity 
peaks with a momentum transfer along the diagonal but 
also the ones with a momentum transfer along the bond 
directions of a square lattice. 





I- J M II-O.J y II 




FIG. 3: The imaginary and real parts of the single particle 
and anomalous Green's function, A(k;E), B(k;E), J(k;E), 
and K{k; E) [as defined by Eqs. ^^-^^], at the energy E = 
— (Aq -|- fio) for tlie electronic coupling to the Big phonon 
(first row), spin resonance mode (second row), and the case 
of no mode coupling. 



The different Fourier spectrum patterns come from the 
detailed renormalization of electronic structure by the 
coupling to these modes. To better understand these 
patterns, we turn to a detailed analytical form of the 
Fourier transform. By putting aside the external form 
factor associated with a specific configuration of weak 
disorder, the Fourier spectrum is determined uniquely 
by the electronic single particle Green's function and is 



found to be: 

5p{ci-E) 
{[A{\^-E)B{\^ + q:,E) 
-[J(k; E)K{\^ + c^, E) ^ 

where 



"0 



-B(k; £;)A(k 
/^(k; i;)J(k4 



q;i?)]},(28) 



A(k; E) ^ Im[gn (k; E + ^7)] , (29) 

TT 

Bi'k-E) - Re[gn(k;S + ^7)] , (30) 
J(k;£;) = --Im[5l2(k;S + ^7)] , (31) 

TT 

K{}i;E) = Re[^l2(k;^; + ^7)] . (32) 

Here as shown in Eqs. (|28|l through 132|1 . the Fourier 
spectrum is determined by the convolution of the imag- 
inarjiS and the real parts of the single-particle {Gii) 
and anomalous (^^12) Green's function in the supercon- 
ducting state. The stronger intensity in 5p{c{] E) (also 
^^'(q;^^)) will be contributed from the wave vector q, 
which connects the simultaneously largest intensity in 
A{k; E) and B(k; E) maps, and in J(k; E) and K{k.] E) 
maps. To be illustrative, we present in Fig.|31those maps 
for the electronic coupling to the Big-1 and spin reso- 
nance modes, and also for the case without the mode cou- 
pling at i? = — (Aq -|- JIq)- Notice that the joint intensity 
of J(k; E)K{k'] E) are smaller by an overall factor of 10 
than that of A{\^\E)B{\^'■,E) and the Fourier spectrum 
is mostly determined by the latter product. For the case 
of no mode coupling, the strongest weight in these quan- 
tities are located at the Al points of the first Brillouin 
zone, which leads to the strongest intensity of the Fourier 
spectrum with the momentum transfer q = (tt, tt) and 
the zero momentum transfer q = 0. If the electrons are 
coupled to the Big-1 phonon mode, besides the strongest 
intensity on the closed ridges in A(k; E) and i?(k; E) [red 
area], there are also moderately strong intensity on the 
two split beams around the M points in these two maps 
[bright green in the former and dark blue in the latter]. 
The intensity A(k; E)B{k';E) connected by q = k' - k 
with k and k' located at the ends of these beams becomes 
stronger. These wave vectors are just those in the Big- 
I panel of Fig.|21at which the Fourier spectrum exhibits 
peaks [green spots]. However, if the electrons are coupled 
to the spin resonance mode, no such beams exist, which 
explains the lack of peaks in Fourier spectrum [see the 
(tt, tt) panel in Fig. [2]. 

Experimentally, only the peaks in the Fourier spectrum 
with the momentum transfer along the bond direction 
has been observed so iarm^ The scenario that the elec- 
trons are coupled to the Big phonon mode and the cou- 
pling matrix is highly anisotropic (i.e., Big-l mode) bears 
the closest resemblance to the experimental observation, 
though not in a full agreement. In Fig.01 we present the 
energy evolution of the spectrum pattern for the elec- 
tronic coupling to the Big-l mode. It shows that the 
characteristic momentum transfer wave vector decreases 
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FIG. 4: The Fourier spectrum of the derivative of the LDOS 
is shown at the various values of the energy for the case of the 
electronic coupling to the Big phonon mode. Here the energy 
has been measured by scaling Aq = 30 meV. 

slightly with the increased energy. However, the intensity 
at these wave vectors decreases rapidly when the energy 
moves away from the action point — (Aq + ^o)- This is 
also not inconsistent with the experiment. 

IV. CONCLUDING REMARKS 

There is considerable evidence that numerous ARPES 
and tunneling data can be interpreted in terms of the 
electronic coupling to bosonic modes with energy about 
40 meV. Possible candidates for this mode are the (tt, tt) 
spin resonance mode and various phonon modes. The 
planar junction tunneling has provided an accurate mea- 
surement of the energy scale of this mode. However, since 
the structure in the tunneling spectra due to the elec- 
tronic coupling to Big phonon mode or to the (tt, tt) spin 
resonance mode does not have much difference, it is very 
difficult to determine the nature of the mode. 

In this paper, we have analyzed the Fourier spectra of 
the energy derivative LDOS to investigate the momen- 
tum transfer structure arising from the electronic cou- 
pling to these different modes. This quantity nicely com- 
plements ARPES for the understanding of the electronic 
responses to these modes. It can now be obtained from 
the FT-STM experiments with elevated spatial resolu- 
tion. In general, we found that the detailed momentum 
dependence of the coupling matrix element strongly in- 
fluences the electronic properties. In particular, wc have 
shown that if the Big or breathing modes are coupled 
to the electrons with only a q dependence, the spectrum 
displays the peak structure along the diagonals with a 
small momentum transfer while the (tt, tt) mode coupling 
does not produce much weight at the small momentum 
transfer region. On the other hand, if the electrons are 
coupled to the Big phonon mode with a matrix element 



that depends not only on q but also on k, the peak struc- 
ture with a small momentum transfer can also appear 
along the Cu-0 bond directions of the Cu02 plane. Re- 
cent FT-STM experiments have indeed found a peak at 
a bond-directed momentum transfer. 

Our calculations also show the structures at large mo- 
mentum transfer for all cases of the electron-collective 
mode coupling. In other words, both types (electronic 
and phononic) of mode couplings produce structures near 
(tt, tt), but only a coupling to the phonon modes yields 
additional structures at the small momentum transfers. 
No peak structures near (tt, tt) have been observed in the 
experiments. The situation is somewhat similar to the 
elastic scattering case (i.e., in the absence of the collec- 
tive mode coupling), where similar structures near (tt, tt) 
also appear in the theoretical spectra£^i£^ but are not ob- 
served experimentally. It is likely that the lack of struc- 
tures at large momentum-transfer in the clastic and in- 
elastic experiments has a common origin. One possibility 
has to do with strong inhomogeneities, which may give 
rise to a dominant forward scattering and make only the 
structures at small momentum transfers observable. This 
amounts to the following form factor for the weak disor- 
der configuration: 



l + r,[sin2(q,/2)+sin2(gy/2)] ' 

where the parameter Tc controls the range of the for- 
ward scattering in the q space. The overall modula- 
tion in the Fourier spectrum of the derivative LDOS 
Sp'{q,uj) (X F{q) will be confined to small momenta if 
F(q) is. 

Finally, several remarks are in order: (i) The electron- 
collective mode coupling we have considered preserves 
the translational symmetry. Fourier spectra were stud- 
ied after introducing clastic impurity scattering with 
weak scattering potential. Alternatively, local vibra- 
tional mode scattering will not only provide an inelastic 
scattering channel but also will break the translational 
symmetry at the beginning»iS. To fully understand the 
possible FT-STM experiments, it would be instructive 
to consider the electronic coupling to the distributed lo- 
cal vibrational modes. We leave this problem for sepa- 
rate analysis, (ii) There has also be increased interest 
in the quasiparticle scattering from the ti impurities^ 
It would be very helpful to study the the FT spectrum 
through the ti scatterers. We delay this investigation to 
a future publication. 
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